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ABSTRACT 

Recent observations show that a large number of Lyman-alpha emitters (LAEs) 
at high redshift z > 3 have unusually high Lya equivalent widths (EW > 400 A). 
However, the origin of these high EWs is an open question. Here, we investigate the 
impacts of photon trapping on the Lya EW and other properties by tracking the Lya 
radiative transfer in spherical galactic clouds. We find that the delayed escape of the 
Lya photons can change the Lya properties significantly. During the transition phase 
from optically thick to optically thin where the Lya photons can escape simultane- 
ously, the EW can be boosted to ~ 1000 A, the Lya luminosity can be increased by 
a factor of a few, and the line profile can be significantly broadened. The boost factor 
appears to depend on the galaxy properties such as mass and star formation rate and 
timcscalc, therefore future investigation combing 3D Lya RT calculations with cos- 
mological simulations of galaxy formation and evolution is needed to fully understand 
the Lya properties of early star-forming galaxies. 

Key words: radiative transfer - galaxies: evolution - galaxies: formation - galaxies: 
high-redshift 



1 INTRODUCTION 

The hydrogen Lya line has played an important role in 
detecting distant galaxies. To date, a large population of 
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mined by the intrinsic spectral energy distribution (SED) of 
the stars. It has been suggested that in the case of a Salpeter 
initial mass fu nction (IMF), the upper limit is ~ 400 A (e.g., 
ISchaererl [20031 ). 

However, a large number of LAEs at z > 3 have 
large EWs highe r than 400 



(e.& 

lOuchi et ahl 120081 '). More recently, Kashikawa et all (|2012D 



Gronwall et all 12007 



portant properties of Lya is the equivalent width (EW), 
which measures the ratio of Lya flux to that of the UV 
continuum, and therefore provides useful information of the 
emission mechanisms and photon escape processes in galax- 
ies. 

Unlike UV continuum, Lya photons typically experi- 
ence many scattering process in galaxies due to its large 
scattering cross section. These scatterings may prolong the 
path length and increase the probability of dust absorption, 
leading to reduction in the emergent luminosity. Hence, the 
EW can change significantly due to the difference in the 
radiative transfer (RT) process between Lya and UV con- 
tinuum. If the escape fraction of Lya is the same as that 
of the UV continuum, the upper limit of the EW is deter- 
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detected a LAE with extremely large EW of ~ 900 A at 
z = 6.5. The origin of such high EWs is largely unknown. 

It has been proposed that the large EWs may be 
produced by a top-heavy IMF of PopII I stars (e.g., 
IGronwall etail 120071 : iKashikawa et all l2012h . However, so 
far there is n o evidence of a top- heavy IMF in high-redshift 
LAEs (e.g., iNagao et all |200S| ), or a PopIII-dominated 
galaxy. Another model used a Salpeter IMF but invoked 
clum py interstellar medium (|Neufeldl Il99ll ; lHansen fc Obi 
I2006D ." In this model, UV continuum is assumed to be ab- 
sorbed by dust passing through the dense clouds, but Lya 
photons do not experience the dust absorption because they 
cannot enter the clouds due to surface scattering. However, 
this is an idealized, optically thin case for the Lya photons. 
In fact, most of the stars form in high density clouds, they 
ionize the gas and create Hu regions which emit Lya pho- 
tons. These photons then travel in the dense cl ouds and may 
experience strong dust extinction. Recently, lYaiima et al.l 
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(|2012bh studied the escape of Lya and continuum photons 
from star-forming galaxies by combining a cosmological sim- 
ulation with three-dimensional radiative transfer calcula- 
tions, and found no correlation between the EW and clumpi- 
ness of the ISM. Fur thermore, using the same simulations, 
lYaiima et al.l (|2012d ) found that collisional excitation from 
accretion of cold hydrogen gas by a galaxy at high redshift 
(z > 6) significantly enhance the Lya emission and boost 
the EW. 

It is not clear, however, from our previous studies 
whether or not cold accretion is the only dominant mecha- 
nism to produce a large EW. In many RT calculations, all 
Lya photons were assumed to escape instantaneously with 
the UV continuum. However, since the ionization of hydro- 
gen evolves with time, this time-dependent propagation of 
the ionizatio n front may cause a d elay of the escape of the 
Lya photons (|Pierleoni et al.f2009l ) , which in turn may affect 
the Lya properties. However, the time evolution of ioniza- 
tion has not been considered in most Lya RT calculations 
before. 

In this work, we consider the pure effects of photon 
trapping on the Lya properties. We perform full Lya RT, 
which includes time-evolving ionization structure, in spher- 
ical galactic clouds, and track the resulting Lya EW, lumi- 
nosity, and line profiles. The paper is organized as follows: 
in § 2 we describe our galaxy models and RT methods, in § 3 
we present the results of the Lya properties, we discuss our 
model assumptions and parameters in § 4, and summarize 
in § 5. 



2 METHODOLOGY 
2.1 Galaxy Models 

Our galaxy model consists of one-dimensional spherical 
shells of hydrogen gas in a spherical N avarro-Frenk- White 
dark matter halo dNavarro et al.l Il997f ) , similar to that of 
iDiikstra fc Loebl (|2009h . The dark matter halo has a con- 
centration parameter c = 3.8 (|Gao et a l. 2008) at th e viria l 
temperature of the halo as derived bv lMakino et al.l (|l998T ) 
at z v ir = 3. The gas is assumed to be in hydrostatic equilib- 
rium and iso thermal, with a t e mpera ture of T = 10 4 K, as 
suggested by IDiikstra fc Loebl (|2009h . A total of 100 shells 
are constructed for each galaxy. 

In addition, we assume star formation takes place at 
the center of the g alaxy, with the rat e SFR proportional to 
the halo mass fe.g. jTrac fe Cer]|2007r i. and decays exponen- 
tially, 

SFR ^^(l^e) exp (-i) W 

where /sfr and tsfr are amplitude factor of peak SFR and 
star formation timescale respectively. In this work, the /sfr 
and tsfr are free parameters with a range /sfr = 0.05 — 5.0 
and tsfr = 5 x 10 6 - 10 8 yr. 

Along with the star formation history, newly formed 
stars are added to the galaxy center. We estimate the num- 
ber of ionizing photons from stacked SEDs. The SED is cal- 
culated by summing that star clust ers are created in e ach 
time step. We use PEGASE ver.2.0 (|Fioc fc Roccalll997l ) to 
calculate the SEDs and assume a metallicity of Z = 0.05Zq. 



In each time step, we estimate the ionization structure from 
the balance between recombination and photo-ionization, 

f Ri 12 

iVion = / 47rr' 2 aBnHii {r)n e {r)dr (2) 
Jo 

where Ni on , Ri, qb are the emissivity of ionizing photons, 
the radius of ionized region and the recombination rate to 
all excitation state, respectively. The gas at r ^ Ri is com- 
pletely ionized region, while at r > Ri it is in neutral state. 



2.2 Radiative Transfer 

The RT calculations are performed using the 3D 
Monte Carlo RT code, All-wavelength Radiativ e Trans- 
fer with Adaptive Refin ement Tree (ART 2 ) (|Li et al.l 
120081 : lYaiima et all l2012al ). The multi-w avelength contin- 
uum emission of ART 2 is described in iLi et ail (|200Sf ). 
w hile the implementati on of Lya line transfer is detailed 
in lYaiima et all (|2012aT ) . Here we focus on the RT of Lya 
photons in the ionization structure of the galactic shells and 
briefly outline the process. 

In the Lya RT, it is important to calculate the fre- 
quency change due to the scattering process. In the scatter- 
ing process, the final frequency in the laboratory frame is 
then: 

^° Ut - VQ = V™ - Up _ \4 ■ fc in V a • fcout 

Aj/ d A^ d V a Vth 1 ' 

where, v ln and v° ut are incoming and outgoing frequency in 
the rest frame of scattering medium, respectively, the line 
center frequency Vo = 2.466 x 10 15 Hz, A^d = (Vth/c)fo 
corresponds to the Doppler frequency width and Vth is the 
velocity dispersion of the Maxwellian distribution describ- 
ing the thermal motions, i.e., Vth = (2kBT/mn) 1 ^ 2 , V a is 
the atom velocity, k- m and fc out are incoming and outgoing 
propagation direction respectively. 

Unlike continuum photons, the Lya photons can expe- 
rience many scattering processes, which lead to long path 
length and long scattering timescale before the photons es- 
cape. As a result, the ionization structure may change while 
the photons are traveling. Here, we consider the change of 
ionization structure with time in the Lya RT process. The 
optical depth is estimated to be dr — a^nm(i, j)ds, where 
i and j are indexes of the shell position and the time step, 
respectively. In each scattering, if the traveling of the Lya 
moves to the next time step, the ionization structure of that 
time step is used to estimate the optical depth. 

The intrinsic Lya luminosity is estimated as: 

L a = 0.68hv a N ion {l - f™) (4) 

where v a and /o°c are the Lya frequency and the escape 
fraction of ionizing photons, respectively. The start point 
of the photon package is randomly chosen with the weight 
by the Lya emissivity of each spherical shell, oc 7r(rf — 
rf_ 1 )nHii^|, where i is the shell index. 

A total of 128 photon packets are used at each time 
step. Each packet has a energy which is the L a at the time 
step divided by the packet number. The initial frequency is 
randomly chosen from a Gaussian distribution. 
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Figure 1. Time evolution of the SFR (upper panel), size of the 
ionization front normalized by the virial radius (middle panel), 
and the escape fraction of ionizing photons (lower panel), for a 
galaxy with mass Mtot = 10 12 Mq , and fiducial parameters star 
formation timescale TgpR = 5 X 10 7 yr, and star formation am- 
plitude factor /sfr = 0.5. 



3 RESULTS 

3.1 Ionization History 

Figure [1] shows the evolution of the star formation rate, 
the size of the ionized region, and the escape fraction of 
ionizing photons with time for a galaxy with mass 

A/tot = 10 12 Mq, a star formation timescale tsfr = 5 x 10 7 
yr, and a star formation amplitude factor /sfr = 0.5. The 
SFR has an initial value of 50 Mq yr -1 , but exponentially 
decreases with time. Since ionizing photons mainly come 
from massive, young stars, they are strongly correlated with 
SFR. As a result, in the early phase (t < 6 x 10 7 yr), the 
ionizing front reaches the R v i T . Then, the ionization front 
steeply decreases with SFR, and becomes ~ 0.1i? v ir at ~ 10 8 

yr- 

In the current model of ionization, the / osc of ionizing 
photons becomes > when all gas in the galaxy is ionized. 
Therefore, the / esc is nonzero at t < 6 x 10 7 yr. It starts from 
~ 0.6, then decreases with SFR. When all gas is ionized, the 
/ csc is estimated by (1 - N ICC}to t)/ N ion , where iVrcc.tot is the 
number of total recombination hydrogen. Since the iVrec.tot is 
time-independent, therefore the / osc is correlated with SFR. 

3.2 Lya Equivalent Width 

We estimate the EW in rest frame from EW = L a /LY V , 
where is luminosity density of UV continuum, here 

we use the luminosity density at 1400 A. Figure [2] shows 
the resulting time evolution of EWs with the fiducial values 
/sfr = 0.5 and tsfr = 5x 10 7 yr, for 4 different galaxies in 
the mass range of 10 10 ~ 13 Mq. 

At first, the ionized region is confined in the galaxy, then 
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Figure 2. Time evolution of the Lya EW for 4 different galaxies 
with mass of 10 10 Mq (upper left panel), 10 11 Mq (upper right 
panel), 10 12 Mq (lower left panel), 10 13 Mq (lower right panel), 
respectively. Fiducial parameters /sfr = 0.5 and tsfr = 5 x 
10 7 yr are used for all galaxies. The red solid line is the resulting 
EW from our model which takes into account the trapping of Lya 
photons due to the time evolution of ionization structure in the 
galaxy, while the black dashed line is the predicted EW assuming 
that all Lya photons escape instantaneously. A Salpeter IMF is 
assumed in these calculations. 



it grows with increasing stellar mass. In this phase, the Lya 
photons are also confined in the galaxy, hence the EW is very 
small. However, when the ionization front reaches the virial 
radius, the trapped Lya photons can escape simultaneously. 
This sudden release of Lya photons causes very high EWs, 
which peaks ~ 400 — 1000 A, much higher than the upper 
limit predicted by the model in which all Lya photons escape 
instantaneously. 

Moreover, the peak of EW decreases with increasing 
halo mass. In small galaxies, the potential well is shallow, 
the SFR is low, and the virial radius R v [ r is small, therefore it 
is easy for the ionization to propagate and reach i? v ir quickly. 
Meanwhile, low SFR produces low UV continuum flux, so 
the resulting Lya EW is high. In contrast, more massive 
galaxies have larger virial radius, so the time it takes for the 
ionization front to reach R V1T is longer, and the higher SFR 
produces higher UV continuum flux, which results in smaller 
Lya EW. In addition, in the larger radius, the traveling time 
from last scattering to escape can change depending on the 
position and angle, resulting in the extended and smaller 
peak of EW. 

Once the galaxy is ionized, the ionizing photons can es- 
cape, and the escape fraction is fixed. From Equations (2) 
and (4), the intrinsic Lya luminosity then becomes constant, 
a.B.n\(r)dr . On the other hand, the 
UV continuum can continue to increase with the formation 
of new stars. Thus, upon completion of ionization, the intrin- 
sic Lya luminosity decreases relatively to UV continuum, 
resulting in a decrease of Lya EW. 
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After that, however, the ionized region is confined 
within _R v ir again due to decreased SFR. Most of Lya pho- 
tons are trapped within the galaxy by the neutral hydro- 
gen shell. Some of the trapped photons can escape through 
scattering and random processes, leading to a mild increase 
of the EW. However, the UV continuum slowly catches up 
again with new star formation, bringing down the EW. 

The interplay between Lya and UV continuum photons 
in the RT processes results in a unique saddle-like shape of 
the EW as a function of time. To better understand the 
oscillatory behavior of the EW, we illustrate in Figure [3] a 
schematic cartoon of the time evolution of the ionization 
structure and the escape of Lya photons. There are four 
major phases in the evolution which affect the Lya EW 
significantly: 

(a) : Early phase when the ionized region is confined in the 
galaxy due to low SFR and stellar mass. Most of the Lya 
photons are trapped in the galaxy. The resulting Lya EW 
is very small. 

(b) : As more stars form, the flood of ionizing photons from 
young, massive stars quickly ionizes the galaxy. The previ- 
ously trapped Lya photons as well as the new ones created 
via recombination of the ionizing photons in the HII regions 
escape instantaneously, causing a sudden and strong boost 
of the EW. 

(c) : With the decline of SFR due to ionization, the number 
of ionizing photons is reduced, and the ionized region is con- 
fined in the galaxy again. The cooled and neutral hydrogen 
can trap the Lya photons, leading to a decrease of the EW. 

(d) : As the stars become older and dimmer, the UV con- 
tinuum flux decreases. Meanwhile, some of the trapped Lya 
photons escape after some scatterings. Hence, the Lya lu- 
minosity increases relatively to the UV continuum, bringing 
a rise of the EW again. 

3.3 Lya Luminosity 

Owing to the time-dependence of the ionization structure, 
the emergent Lya luminosity (L a ) from the galaxy evolves 
with time as well. Figure 0] shows the corresponding L a from 
the modeled galaxies as a function of time. 

During Phase (a), the Lya photons are trapped in the 
galaxy, the L a approaches to zero in the modeled galax- 
ies. As the galaxy gradually becomes optically thin due to 
ionization, Lya photons can escape, steadily increasing the 
luminosity. Upon Phase (b) when the galaxy is completely 
ionized, all the Lya photons, including both the previously 
trapped and the newly created ones, escape instantaneously, 
resulting in a drastic, strong peak of the L a . After that, the 
ionized region is confined in the galaxy, the Lya photons 
are trapped again as in Phase (c) and (d). But since the 
Lya photons escape mainly through scattering, the fraction 
remains nearly constant during these phases, leading to con- 
stant luminosity. At the end of the star formation timescale, 
no more new stars form, the number of ionizing and Lya 
photons declines sharply, so the galaxy becomes optically 
thick, and most of the Lya photons become trapped, result- 
ing in a rapid drop of the luminosity. 

The evolution of the Lya luminosity does not show the 
characteristic "saddle" shape as the EW. This is because the 
Lya luminosity depends only on the ionization structure and 




log t (yr) 



Figure 4. Evolution of emergent Lya luminosity from the four 
modeled galaxies as indicated by the colors (cyan : Aftot = 
10 10 Af Q , blue : Aftot = 10 11 Mq, red : Aftot = 10 12 Af and 
green : Af to t = 10 13 Af Q ). 



the escaping process, it does not depend strongly on the UV 
continuum as does the EW. 

In contrast to the EW, which in general decreases with 
galaxy mass, the Lya luminosity increases with galaxy mass. 
This is because both the SFR and the total recombination 
number are proportional to the galaxy mass, and as a result, 
the L a increases linearly with the galaxy mass. 

3.4 Lya Line Profiles 

Figure \5\ shows the emergent line profiles from the model 
which takes into account the delayed escape of Lya photons 
due to the time-depdendence of the ionization structure in 
the galaxy, in comparison with those without photon trap- 
ping. 

At t = 6.5 x 10 5 yr, which corresponds to Phase (a) 
in Figure [3j the delayed-escaping model predicts that the 
Lya photons are trapped in the galaxy, there is no emergent 
line emission. In contrast, the instantaneous-escaping model 
predicts a double-peak profile caused by scattering of the 
neutral shells. At Phase (b) when the galaxy is completely 
ionized (t = 1.2 x 10 6 yr), the galaxy becomes optically thin. 
The line profile without the time dependence shows a narrow 
profile peaking around x = 0, which is the intrinsic profile, 
while the line from the trapping model shows not only the 
central peak, but also a significantly extended distribution 
in both the blue and red wings. The broadening is caused by 
the scatterings of the trapped photons before they escape. 
Then, in Phase (c) (t = 9.0 x 10 7 yr), the galaxy becomes 
optically thick again, the line profile thus shows double peak 
owing to scattering, and the luminosity from the trapping 
model is lower than that from the no-trapping model as ex- 
pected. In Phase (d) (t ~ 1.3 x 10 s yr), with the decline of 
SFR, the size of the ionized region becomes small. As a re- 
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(a) (b) (c) (d) 




Figure 3. Schematic view of the interplay between ionization and Lya escaping as a function of time, (a): the ionized region is confined 
in the galaxy, the Lya photons are trapped, (b): the galaxy is completely ionized, the accumulated Lya photons escape instantaneously, 
(c): the ionized region is confined again, and the Lya photons are trapped again, (d): the stars become older and dimmer, some Lya 
photons escape through scattering. 
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Figure 5. The Lya line profile of the Mtot = 10 12 Mq galaxy 
at four different times: t = 6.5 X 10 5 , 1.2 X 10 6 , 9.0 X 10 7 and 1.3 X 
10 s yr, respectively. Fiducial parameters star formation timcscale 
tsfr = 5 X 10 7 yr, and star formation amplitude factor /sfr = 
0.5 are used. The red solid line is the line profile from our model 
which takes into account the trapping of Lya photons due to 
the time evolution of ionization structure in the galaxy, while 
the black dashed line is the predicted profile assuming that all 
Lya photons escape instantaneously. The x is frequency shift in 
Doppler units, x = (y — ^o)/Ai^rj. 



suit, the optical depth from the center to the 7? v ir increases, 
the distance between the peaks in the profile thus increases, 
resulting in two broad double peaks. 



4 DISCUSSIONS 

In our model, there are two free parameters, the amplitude 
factor of peak SFR, /sfr, and star formation timescale, 
tsfr. In the previous sections, we have presented results 
using the fiducial numbers /sfr = 0.5 and tsfr = 5 x 10 7 
yr. Here we explore the dependence of our results on these 
two parameters. 




log t (yr) log t (yr) 

Figure 6. The EW with different amplitude factor of SFR as a 
function of evolution time. The total mass of galaxy is 10 12 Mq 
and tsfr = 5x 10 7 yr. Black dash lines are the EW in the assump- 
tion that all Lya photons instantaneously escape. Red solid lines 
are the EW considering the time evolution of ionization structure 
while Lya photons travel in the galaxy. 



Figure |B] shows the EW with difference /sfr for the 
galaxy with Mtot = 10 12 M & and t SF r = 5 x 10 7 yr. For 
/sfr = 0.05, the ionized region is always confined in the 
galaxy owing to low SFR. Hence, there is no EW peak even 
at the Phase (b) in figure [3] The EW simply decreases with 
SFR as the mean stellar age becomes old. For /sfr ~ 0.5 — 
1.0, there is a EW peak at Phase (b). With higher /sfr, 
however, both the SFR and the mass of young stars are 
higher, leading to an earlier completion of ionization and 
higher UV continuum flux, which results in smaller EW. 

Figure [7] shows the EW with different tsfr for the 
galaxy with M to t = 10 12 Mq and /sfr = 0.5. All plots 
show significantly boosted EW peak at Phase (b), and sim- 
ilar "saddle" shape of the evolution of EW. The difference 
among these models lies in the duration of Phase (c). With 
shorter tsfr, the the occurrence of Phase (b) is more sig- 
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Figure 7. The EW with different time scale of star formation as 
a function of evolution time. The total mass of galaxy is 10 12 Mq 
and /sfr = 0.5. Black dash lines are the EW in the assumption 
that all Lyo photons instantaneously escape. Red solid lines are 
the EW considering the time evolution of ionization structure 
while Lyo photons travel in the galaxy. 



nificantly delayed, resulting in larger difference in the EW 
evolution between the trapping and non-trapping models. 

These results suggest that the EW depend sensitively 
on the galaxy properties such as star formation rate and 
timescale. In order to fully address the origin and nature 
of the EW distribution in LAEs, it is desired to combine 
3D Lya RT calculations with full cosmological simulations 
which include more realistic treatments of star formation 
and feedback processes, which we plan to do in future work. 



• Phase (d): The ionized region becomes smaller, as the 
stars become older and dimmer. The escape of the trapped 
Lyo photons is dominated by the stochastic scattering pro- 
cess. There is a small enhancement of the EW due to the 
decrease of the UV continuum flux, but the Lya luminosity 
remains constant until the galaxy becomes strongly optically 
thick and trap most of the photons. 

Our model suggests that the interplay between time- 
dependent ionization and Lya scattering impacts signifi- 
cantly the Lya EW, the emergent luminosity, and the line 
profile. During the transition phase from optically thick to 
optically thin where the Lya photons can escape simultane- 
ously, the EW can be boosted to ~ 1000 A, the Lya luminos- 
ity can be increased by a factor of a few, and the line profile 
can be significantly broadened. These results may be able 
to explain the unusually large EW in some observed LAEs. 
However, the boost factor depends on the galaxy proper- 
ties such as mass and star formation rate and timescale, 
future investigation combing 3D Lya RT calculations with 
cosmological simulations of galaxy formation and evolution 
is needed to fully address this issue. 
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5 SUMMARY 

We have performed a set of Lya RT calculations on idealized 
spherical galaxies to study the effects of photon trapping on 
the Lya properties in star-forming galaxies, by taking into 
account the time dependence of the ionization process. We 
have identify four major phases in the evolution of the galaxy 
in our model: 

• Phase (a): The ionized region is confined in the galaxy 
due to low SFR and stellar mass. Most of the Lya photons 
are trapped in the galaxy, no emergent Lya emission is ob- 
served. 

• Phase (b): The galaxy is completely ionized. The pre- 
viously trapped Lya photons as well as the new ones cre- 
ated via recombination escape instantaneously, resulting in 
a dramatic and strong boost of the EW and luminosity, and 
a broadened line profile. 

• Phase (c): The ionized region is confined in the galaxy 
again due to decrease of SFR and ionizing photons. The Lya 
photons are trapped again, leading to a decrease of the EW 
and luminosity. 
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